Lab9
In [1]:
Copied!
import gtlparser
import gtlparser
Question 1: Creating an Interactive Map¶
In [2]:
Copied!
map1 = gtlparser.Map()
map1.add_basemap_gui(position="topright")
map1.add_layer_control()
map1.add_search_control(position="topleft")
map1
map1 = gtlparser.Map()
map1.add_basemap_gui(position="topright")
map1.add_layer_control()
map1.add_search_control(position="topleft")
map1
Out[2]:
Question 2: Adding Map Legends¶
In [3]:
Copied!
import leafmap
map2 = leafmap.Map()
wms_url = "https://services.terrascope.be/wms/v2?"
layer_name = "WORLDCOVER_2021_MAP"
map2.add_wms_layer(
url=wms_url,
layers=layer_name,
name="ESA World Cover 2021",
format="image/png",
transparent=True,
shown=True,
version="1.3.0",
)
map2.add_legend(builtin_legend="ESA_WorldCover")
map2
import leafmap
map2 = leafmap.Map()
wms_url = "https://services.terrascope.be/wms/v2?"
layer_name = "WORLDCOVER_2021_MAP"
map2.add_wms_layer(
url=wms_url,
layers=layer_name,
name="ESA World Cover 2021",
format="image/png",
transparent=True,
shown=True,
version="1.3.0",
)
map2.add_legend(builtin_legend="ESA_WorldCover")
map2
Out[3]:
Question 3: Creating Marker Clusters¶
In [4]:
Copied!
import leafmap
import geopandas as gpd
import pandas as pd # geopandas relies on pandas
geojson_url = "https://github.com/opengeos/datasets/releases/download/places/wa_building_centroids.geojson"
gdf = gpd.read_file(geojson_url)
gdf
gdf["longitude"] = gdf.geometry.x
gdf["latitude"] = gdf.geometry.y
map3_1 = leafmap.Map()
map3_1.add_marker_cluster(gdf, x="longitude", y="latitude", layer_name="World cities")
map3_1
import leafmap
import geopandas as gpd
import pandas as pd # geopandas relies on pandas
geojson_url = "https://github.com/opengeos/datasets/releases/download/places/wa_building_centroids.geojson"
gdf = gpd.read_file(geojson_url)
gdf
gdf["longitude"] = gdf.geometry.x
gdf["latitude"] = gdf.geometry.y
map3_1 = leafmap.Map()
map3_1.add_marker_cluster(gdf, x="longitude", y="latitude", layer_name="World cities")
map3_1
Out[4]:
In [5]:
Copied!
map3_2 = leafmap.Map()
map3_2.add_circle_markers_from_xy(
gdf,
x="longitude",
y="latitude",
radius=5,
color="red",
fill_color="yellow",
fill_opacity=0.8,
marker_cluster=True,
)
map3_2
map3_2 = leafmap.Map()
map3_2.add_circle_markers_from_xy(
gdf,
x="longitude",
y="latitude",
radius=5,
color="red",
fill_color="yellow",
fill_opacity=0.8,
marker_cluster=True,
)
map3_2
Out[5]:
Question 4: Visualizing Vector Data¶
In [6]:
Copied!
geojson_url = "https://github.com/opengeos/datasets/releases/download/places/wa_overture_buildings.geojson"
map4_1 = leafmap.Map()
geojson_style = {
"color": "red",
"weight": 2,
"fillOpacity": 0,
}
map4_1.add_geojson(geojson_url, style=geojson_style, layer_name="Washington Buildings")
map4_1
geojson_url = "https://github.com/opengeos/datasets/releases/download/places/wa_overture_buildings.geojson"
map4_1 = leafmap.Map()
geojson_style = {
"color": "red",
"weight": 2,
"fillOpacity": 0,
}
map4_1.add_geojson(geojson_url, style=geojson_style, layer_name="Washington Buildings")
map4_1
Out[6]:
In [7]:
Copied!
geojson_url = "https://github.com/opengeos/datasets/releases/download/places/las_vegas_roads.geojson"
map4_2 = leafmap.Map()
geojson_style = {
"color": "red",
"weight": 2,
}
map4_2.add_geojson(geojson_url, style=geojson_style, layer_name="Las Vegas Roads")
map4_2
geojson_url = "https://github.com/opengeos/datasets/releases/download/places/las_vegas_roads.geojson"
map4_2 = leafmap.Map()
geojson_style = {
"color": "red",
"weight": 2,
}
map4_2.add_geojson(geojson_url, style=geojson_style, layer_name="Las Vegas Roads")
map4_2
Out[7]:
In [8]:
Copied!
map4_3_data = (
"https://github.com/opengeos/datasets/releases/download/us/us_counties.geojson"
)
map4_3 = leafmap.Map(center=(40, -98), zoom=4)
map4_3.add_data(
map4_3_data,
layer_name="US County Area (Census)",
column="CENSUSAREA",
scheme="EqualInterval",
cmap="Blues",
legend=True,
legend_title="Census Area (sq mi)",
)
map4_3
map4_3_data = (
"https://github.com/opengeos/datasets/releases/download/us/us_counties.geojson"
)
map4_3 = leafmap.Map(center=(40, -98), zoom=4)
map4_3.add_data(
map4_3_data,
layer_name="US County Area (Census)",
column="CENSUSAREA",
scheme="EqualInterval",
cmap="Blues",
legend=True,
legend_title="Census Area (sq mi)",
)
map4_3
Out[8]:
Question 5: Creating a Split Map¶
In [9]:
Copied!
pre_event_url = (
"https://github.com/opengeos/datasets/releases/download/raster/Libya-2023-07-01.tif"
)
post_event_url = (
"https://github.com/opengeos/datasets/releases/download/raster/Libya-2023-09-13.tif"
)
map5 = leafmap.Map()
leafmap.split_map(left_layer=pre_event_url, right_layer=post_event_url)
map5
pre_event_url = (
"https://github.com/opengeos/datasets/releases/download/raster/Libya-2023-07-01.tif"
)
post_event_url = (
"https://github.com/opengeos/datasets/releases/download/raster/Libya-2023-09-13.tif"
)
map5 = leafmap.Map()
leafmap.split_map(left_layer=pre_event_url, right_layer=post_event_url)
map5
Out[9]: